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Ordinary least square (OLS) and robust least square (RLS) consisting of 
least absolute residual and Bi-square approaches were deployed to obtain the 
fuel consumption characteristic curve and the coefficients of the quadratic 
fuel consumption function for thermal stations in Nigeria. Results were 
compared based on convergence property, root mean square error, R-square 
value, adjusted R-square value, and width interval of coefficients. Valve 
Point loading effects of Egbin and Sapele power stations were used to 
develop the quadratic fuel consumption characteristic curve and function. 
The average difference in width interval for the coefficients a, b, c, d, and e 
of the two stations, after comparing Bi-square and OLS technique, were 
0.02084, 8.5055, 1856.565, 520.8855, and 0.0082, respectively. The R- 
square values obtained from the Bi-square technique were superior to the 
OLS technique with arithmetic differences of 0.6196 and 0.5254 for Egbin 
and Sapele generating stations, respectively. Bi-square technique also 
offered better adjusted R-square value for Egbin and Sapele stations with 
arithmetic differences of 0.622 and 0.5287, respectively. Bi-square technique 
offered smaller root mean square error than the ordinary least square 
technique for both stations. The coefficients obtained from Bi-square 
technique were used to develop the fuel quantity function for the studied 
stations. 
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1. INTRODUCTION 


The operating cost is a vital aspect of electrical energy production and it constitutes the economics of 
the operation. The optimal solution to power flow and economic dispatch problems have become critical and are 
now considered research issues in modern power systems [1]-[3]. These problems are formulated 
mathematically and their solutions are optimize based on fuel cost function subject to various operational 
constraints [4]-[7]. The solution obtained from solving these problems will save generating costs, especially in 


fossil fuel plants. 


The fuel cost function which is commonly represented by a quadratic function is affected by many 
factors vis, the ambient temperature of operation and the rate of aging of the generating units. It is a significant 
issue to accurately estimate the thermal unit fuel cost curve parameters periodically [8]. Therefore, using a 
reliable and powerful estimation technique in the estimation of parameters of fuel cost curve is very imperative. 
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Mathematical methods have become a vital technique for estimating coefficients and data fitting. It has 
become a measure for statistical inference [9]. Till date, only few works have been done on the estimation of 
fuel cost curve parameters. In the past decades, several techniques were employed in solving the complex fuel 
cost problem. Each method has its own advantages and disadvantages. Different mathematical models were 
used to represent the fuel cost function [8]. The basic principle of the OLS method is to estimate the parameter 
vectors, by minimizing the sum of squares of errors in order to find the best function that matches a set of 
data [10]-[12]. Least square techniques aim at reducing the sum of squared errors between the obtained data and 
the actual data and has been used for state estimation in power system [13]-[16]. Taylor’s series expansion 
principle of the least squares error curve fitting technique was used to estimate power system frequency and 
amplitude [17]. Least square technique has found diverse applications in electrical engineering such as in power 
quality assessment [18], determine symmetrical components of power system [19], phasor estimation for power 
quality monitoring [20], to mention but a few. 

Economic dispatch problem is a problem of power system management that can be solved when the 
fuel quantity function of thermal plants is obtained. The ordinary least-square (OLS) method was used to 
develop fuel quantity functions of thermal generating plants [21], [22]. Although, the accuracy of OLS is 
threatened by the presence of outliers in data set [23]. The need to explore the benefits of robust techniques for 
the development of quadratic fuel function of thermal plants becomes key. In order to optimally manage thermal 
generating stations, the coefficients of quadratic fuel quantity function must be uniquely determined. Most 
literature obtained these coefficients using the ordinary least square technique that is sensitive to extreme values. 
This has motivated me to use the robust technique. This study is aimed at developing the quadratic fuel quantity 
function of steam power stations with a focus on Egbin and Sapele power stations both in the Nigeria power 
grid. The objective of this work include the development of the coefficients (a, b, c, d, e) of the fuel quantity 
function while considering the effect of valve point loading and accessing the goodness-of-fit of function 
developed. 


2. FORMULATION OF LEAST SQUARE MODEL FOR FUEL QUANTITY FUNCTION OF 
THERMAL STATION 

Consider the polynomial to degree n for a functional relationship between x and y. 

y = ao +a,x + a3x? + + apx” (1) 

y = Vio a x! (2) 
For a set of N-data points (x;,y;), given a set of points (x1 , Y1),(X2 , Y2),(X3 , Y3)». --, (Xn Yn) the least square 
estimates of ao, a1, a2 ,...,a, can be obtained in a similar way as that of linear regression. The sum of 
squared deviation of the values of y from predicted value is given by [24]. 

S = (y — ao — 41X — a2X? nx")? (3) 
In (3) can be minimized by setting its partial derivatives with respect to (4). 

ao, Ay, Az, ..,A, = 0 (4) 
The following equation is obtained: 

Ly =a +a, Dx +a, yx?...a, Dx (5) 
the cost function for the thermal stations in Nigeria can be described with the (6) [25], [26]. 

Fi (Pei) = Li (aiPGi + BiPoi + ci) (6) 
Where a;, b;, and c; are coefficients of the it” generating station. n is the total number of thermal stations and 
Pg is the power generated by the it” generating station. 


For Valve Point Loading effect consideration which is the case of Egbin and Sapele steam power 
stations, in (6) can be modified thus: 


F,(Pgi) = Lier (aiPG; + b:Poi + ci + 


di sin (e,(Pm" — P:))|) (7) 
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hence, the sum square error is given by (8). 
s = YL — (aP; + bP; + c + |dsin(e(P™™ — PD) |)? (8) 


In order to develop optimal coefficients for the objective fuel cost function of the steam thermal stations, the 
robust least square algorithm (Bi-square algorithm) is used to minimize (8) via the following steps: 

1) Fit the model by weighted least squares; 

2) Compute the adjusted residuals and standardize them. The adjusted residuals are given by (9); 


Tadjusted = TE (9) 
ri are the usual least-squares residuals and z; are leverages that adjust the residuals by reducing the 
weight of high-leverage data points having a large effect on the least-squares fit. The standardized 
adjusted residuals are given by (10); 


Tag: 
u = adjusted (10) 
cs 


C is the tuning constant=4.685 
S is the robust standard deviation 
3) Compute the robust weights as a function of u. The Bi-square weights are given by (11); 


maf GOD hal < i 


0 Ju;| = 1 


4) Perform the next iteration of the fitting procedure by returning to the first step if the computation fails to 
converge at first iteration. 


3. DEVELOPMENT OF FUEL QUANTITY FUNCTION 

The fuel quantity consumed by the stations and their power output are compiled using Microsoft 
Excel Package. This data which constitutes the input-output data for the stations is used to develop the 
input-output characteristic curve using curve fitting tools in MATLAB 2017 model. The coefficients of fuel 
quantity function are obtained from the curve alongside the ‘goodness of fit’. 


3.1. Egbin power station 

The relationship between the quantity of fuel utilized by the station and the power output from the 
generator (input-output relationship) is described by the graph in Figure 1. This relationship curve is obtained 
using the ordinary least square approach (OLS). The goodness of fit is described thus: R-square=0.3651; 
adjusted R-square=0.3626; RMSE=2364. The residual plot can be seen in Figure 2. This plot shows how 
much the outliers differ from the fitted model. 

The relationship between the fuel input and power output for the generating station is obtained using 
a robust least square approach known as Bi-square. The relationship curve is shown in Figure 3. The 
goodness of fit is described thus: R-square=0.9847; adjusted R-square=0.9846; RMSE=367.1. The residual 
plot can be seen in Figure 4. 

After the test carried out on the three approaches (OLS and Bi-square), the Bi-square approach 
produced the best goodness of fit. It had a lower RMSE, a higher R-square and a higher adjusted R-square 
than the OLS approach. Hence, the coefficients together with the width intervals obtained from the Bi-square 
approach are thus; 
a=4.288e-05 (0.0010728) 
b=12.09 (1.64) 
c=11.1 (616.6) 
d=6.447 (150.43) 
e=0.1453 (0.0152) 

The quadratic fuel quantity function for Egbin power generating station is given thus: 


F,3(Py3) = 0.00004288P?, + 12.09P,3 + 11.1 + |6.447sin (0.1453(PI3'" — P,;))| (12a) 
from the curve fitting, 
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pmin = 68.9MW 

hence, (8a) can be modified thus; 
F,3(Py3) = 0.00004288P?, + 12.09P,3 + 11.1 + |6.447sin (0.1453(68.9 — P,3))| (12b) 


these coefficients were computed with 95% confidence bound. 
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Figure 1. Input-output curve for Egbin power generating station (OLS) 
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Figure 2. The residual plot from the input-output characteristics curve of Egbin power station 
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Figure 3. Input-output curve for Egbin power station (Bi-square) 
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Figure 4. The residual plot from the input-output characteristics curve of Egbin power generating station 


3.2. Sapele steam power generating stations 

The relationship between the quantity of fuel utilized by the station and the power output from the 
generator (input-output relationship) is described by the graph in Figure 5. This relationship curve is obtained 
using the ordinary least square approach (OLS). The fuel input which is measured in MMBTU/hr is 
represented by the vertical axis (h14). The power output from the station is represented by the horizontal axis 
(P14). The goodness of fit is described thus: R-square=0.453; adjusted R-square=0.4496; RMSE=296.1. 

The residual plot can be seen in Figure 6. This plot shows how much the outliers differ from the 
fitted model. The relationship between the fuel input and power output for the generating station is obtained 
using a robust least square approach known as Bi-square. The relationship curve is shown in Figure 7. The 
goodness of fit is described thus: R-square=0.9784; adjusted R-square=0.9783; RMSE=58.77; The residual 
plot can be seen in Figure 8. 

After the test carried out on the two approaches (OLS and Bi-square), the Bi-square approach 
produced the best goodness of fit. It had a lower RMSE, a higher R-square and a higher adjusted R-square 
than the OLS approach. Hence, the coefficients as well as the width intervals obtained from the Bi-square 
approach are thus: 
a=0.002837 (0.0089360) 
b=8.989 (2.022) 
c=23.29 (105.27) 
d=18.38 (28.739) 
e=0.5412 (0.01) 

The quadratic fuel quantity function for Sapele Steam power generating station is given thus; 


F,4(Py4) = 0.002837P2, + 8.989P,, + 23.29 + |18.38sin (0.5412(PH" — P,,))| (13a) 
from the curve fitting, P77" = 23.48MW 


hence, (9a) can be modified as (13b). 


Fy4(P,4) = 0.002837P2, + 8.989P,, + 23.29 + |18.38sin (0.5412(23.48 — Py4))| (13b) 
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Figure 5. Input-output curve for Sapele steam power generating station (OLS) 
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Figure 6. The residual plot from the input-output characteristics curve of Sapele steam power station 
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Figure 8. The residual plot from the input-output characteristics curve of Sapele steam power generating 
station 


The width interval of the fuel quantity coefficients for the generating stations is contained in 
Table 1. Width interval is another indicator for goodness of fit. A smaller interval indicates a better fit. From 
Table 1, it can be seen that bi-square offers a narrower width interval. After comparing both techniques for 
Egbin generating station, the differences in width interval (OLS value-Bi-square value) for the various 
coefficients were as follows; a-coefficient was 0.0058102, b-coefficient was 8.863, c-coefficient was 3288.4, 
d-coefficient was 919.87 and e-coefficient was 0.0056. The differences in width interval for the various 
coefficients obtained by comparing both techniques for Sapele steam station were as follows; a-coefficient 
was 0.035864, b-coefficient was 8.148, c-coefficient was 424.73, d-coefficient was 121.901 and e-coefficient 
was 0.0108. After considering the difference in width interval for the coefficients of both stations, average 
values were obtained thus; a-coefficient was 0.02084, b-coefficient was 8.5055, c-coefficient was 1856.565, 
d-coefficient was 520.88550 and e-coefficient was 0.0082. The implication of this result is that the 
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coefficients obtained from the Bi-square technique have a lower degree of uncertainty as compared with the 


results from OLS because they have a narrower width interval. Therefore, Bi-square technique produced a 
better fit than the ordinary least square technique. 


Table 1. Width interval of fuel quantity coefficients for generating stations (hydro-thermal plant) 


S/No. Generating Station A B C D z 
: g Bi-Square OLS _Bi-Square OLS _Bi-Square OLS Bi-Square OLS _Bi-Square OLS 
1 Egbin 0.0010728 0.006883 1.64 10.503 616.6 3905 150.43 1070.3 0.0152 0.0208 


2 Sapele Steam 0.008936 0.0448 2.022 10.17 105.27 530 28.739 150.64 0.01 0.0208 


A summary table showing the goodness of fit parameters for the curve fitting of the thermal stations 
is shown in Table 2. The arithmetic differences between the R-square values obtained from the two 
techniques were 0.6196 and 0.5254 for Egbin and Sapele power stations respectively. The arithmetic 
differences between the adjusted R-square values obtained from the two techniques were 0.622 and 0.5287 
for Egbin and Sapele power stations respectively. The Bi-square technique offered a better ‘goodness of fit’ 
parameters than the ordinary least square approach. 

The R-square and adjusted R-square values obtained from the Bi-square technique were better than 
that of the OLS technique for Egbin and Sapele generating stations because the values for the Bi-square 
technique were closer to 1 than that of the OLS technique. The RMSE values obtained from the Bi-square 
technique were less than that from OLS technique hence, Bi-square technique offered a better fit. 


Table 2. Goodness of fit for the curve fitting of thermal plants (steam power plants) 


; was R-Square Adjusted R-Square RMSE 
penn Generating Station Bi-Square OLS __Bi-Square OLS __Bi-Square __OLS 
1 Egbin 0.9847 0.3651 0.9846 0.3626 367.1 2364 
2 Sapele Steam 0.9784 0.4530 0.9783 0.4496 58.77 296.1 


4. RESULTS 

It was seen that the the width intervals for the various coefficients were narrower with the Bi-square 
technique hence indicating that the degree of uncertainty in obtaining the coefficients was low as compared 
with the OLS technique. It was also observed that the R-square and adjusted R-square values (coefficient of 
determination) obtained by using the Bi-square technique were closer to 1 than that of the OLS technique. 
The RMSE values obtained from the Bi-square technique were less than that of the OLS technique. 


5. CONCLUSION 

This work compared two classes of techniques used for the development of the quadratic fuel 
quantity for the different thermal stations in Nigerian 330 kV power system network. The Bi-square approach 
which is a robust least square technique offered better performance in terms of its goodness of fit when 
compared with the ordinary least square approach. The Bi-square algorithm produced the best fitting when 
compared with the ordinary least square approach. 

In the case of Egbin power station and sapele steam power station, the Bi-square technique produced 
a better width interval for coefficients a, b, c, d, and e. The difference in width interval obtained by 
comparing both techniques for Egbin power station is as follows: The difference in width interval obtained 
by comparing both techniques for Egbin power station produced the following values; 0.0058102, 8.863, 
3288.4, 919.87 and 0.0056 for a, b, c, d, and e coefficients respectively. The difference in width interval 
obtained by comparing both techniques for Sapele steam power station gave a value of 0.035864, 8.148, 
424.73, 121.901, and 0.0108 for a, b, c, d, and e coefficients respectively. These values were obtained based 
on the width interval of fuel quantity coefficients for generating stations (hydro-thermal plant). 

The Bi-square offered better R-square for Egbin and Sapele steam power stations with an arithmetic 
difference of 0.6196 and 0.5254 respectively. Bi-square also offered better adjusted R-square value for Egbin 
and Sapele power stations with a difference of 0.622 and 0.5287. Bi-square also produced smaller root mean 
square error than the ordinary least square technique for both stations. The arithmetic differences were 
obtained using the goodness of fit for the curve fitting of thermal plants (steam power plants). The valve 
point loading effect of Egbin power station and Sapele steam power stations were considered during the 
development of the quadratic fuel quantity characteristics. The quadratic fuel quantity characteristics 
developed can be used as building blocks for carrying out economic dispatch operations. 
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